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摘要 : 本 文 的 目的 研究 磁 隧 减 对 具有 两 成 分 模型 的 脉冲 星 目 转 减 速 的 作用 。 利用 


分 析 方 法 研究 了 两 成 分 模型 的 脉冲 星 在 磁 衰 减 制 劲 力矩 作用 下 两 成 分 的 自转 角 


速度 随时 间 变 化 。 给 出 了 具有 磁 衰 减 两 成 分 模型 的 耦合 方程 组 和 其 解 。 利 用 分 析 


解 对 具有 两 成 分 的 蟹 状 星云 脉冲 星 (PSR0531+21，Crab) 和 船 帆 座 脉冲 星 


(PSR0833-45, Vela) 在 磁 衰 减 作 用 下 做 了 数值 计算 。 数 值 结 果 表 明 : 两 个 脉冲 


星 的 自转 角速度 逐年 减 慢 每 年 为 Crab WEN -0.1710704/), Ti Vela 减速 为 


-0.007 人 ad/)。 最 后 讨论 了 文中 所 得 到 的 结果 并 给 出 在 两 成 分 模型 中 磁 衰 减 在 


在 的 结论 。 
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由 


己 知 仍 状 星云 脉冲 星 (Crab〉 和 船 帆 座 脉冲 星 (Vela〉 除 有 星 震 外 还 具有 两 


成 分 模型 。 外 层 是 带电 的 固体 外 这， 内 层 是 中 子 超 流体 。 两 成 分 在 厅 合 力矩 作用 


下 两 者 以 不 同 角速度 绕 目 转 轴 旋 转 。 和 耦合 力矩 有 引力 辐射 力矩 ， 磁 辐射 力 算 和 磁 


衰减 力矩 中 。 作 者 在 前 文 [4] 研究 了 磁 辐 射 制 动力 矩 对 艇 状 星云 脉冲 星 (Crab) 


的 自转 角速度 减速 的 影响 ,本文 继 前 文 研究 了 磁 辟 减 制 动 力矩 对 稻 状 星云 脉冲 星 


(Crab) 和 有 船 帆 座 脉 冲 星 〈Vela) 的 自转 角速度 的 减速 影响 。 
1. 在 磁 衰 减 制 动 力矩 作用 下 的 两 成 分 模型 的 方程 组 
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作者 在 文 [4] 根 据 文 [1] 给 出 的 两 成 分 模型 在 看 合力 矩 作用 下 的 基本 方程 组 : 


I 
be says). (1) 
dt T 
I, Ne 2G Q,). . (2) 
dt T, 
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万 分 别 是 内 外 两 层 的 转动 惯量 : N 是 作用 在 两 成 分 的 力矩 ，z 是 微观 的 驰 预 时 


间 ; 7 是 宏观 驰 预 时 间 。 它 们 之 间 的 关系 在 文 [4] 中 已 经 给 出 


1=1,4+1, . (3) 


其 中 工 是 两 成 分 的 总 转动 惯量 


中 子 超 流体 丰富 度 Q 
Q=", 1,=01 (4) 
aT a0 . (5) 
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以 下 推出 在 磁 衰 减 制 动 力 窍 作用 下 的 力矩 形式 。 根 据 脉 冲 星 的 磁 偶 极 辐射 模型， 


前 文 [和 利用 
dQ 203 
ae a a 


AP u AWE, u=R BSina ，R HIEP, BIRM. ERRATE 


HAH, a=90, u=R°B. 


正如 前 文 [将 上 面 方程 改 为 角 动 量 SOT RE 


2 
E E per (9) 
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J(t0) 是 to 时 的 角 动 量 。 如 果 取 脉冲 星 目前 (现在) 的 角速度 和 角 动 量 , Bll = 0 时 


的 角 动 量 ， 上 述 积 分 取 上 下 限 为 


六 dJ 2p a 


-= f (10) 
10 F 3 *(CT) 


磁 矩 妈 随 时 间 的 衰减 用 指数 形式 品 ， 则 

(Ê =p e”. wo=RB, . (11) 
<“ 为 磁 衰 减 系数 。 
将 (11) 式 的 4 代入 (10) 式 积 分 后 有 


10D=7Oh-K(e* -jj (12) 


Atty Oo 


式 中 常数 K = ee 


(13) 


(12) 式 和 和 (13) 式 同 前 文 [4] 的 J 和 KK 大 有 不 同 。 
由 (12) 式 给 出 力矩 N 的 形式 : 


N= - osh K(e* -1)[”, 


wee STOK" h- zle” 1)”. (14) 


根据 前 文 [4] 转 动 惯量 1. 和 不 变 情况 下 ，(12) 式 的 角 动 量 


IOn =1.Q+1,0,=7(0),h- Kle -1)). (15) 
将 上 式 改写 为 


0, =h- Ke-1] -0 (16) 


n n 


式 中 J(0), =10(0), FI (4) 式 和 (6) 式 


15e A 


将 其 代入 (16) 式 后 


o, =- -xet -4 - 2a (17) 
Q Q 
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再 将 C14) 中 的 N 和 (17) 中 的 9, 代入 (1 式 , 再 利用 (7) 式 二 = Loe 
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AH (7) 式 


内 部 中 子 超 流体 层 在 人 磁 衰 减 制 动 作用 下 自转 角速度 方程 式 


Wy Oo -2g . (19) 
dt T T 


2. 在 磁 衰 减 制 动力 矩 作用 下 两 成 分 耦合 方程 组 的 分 析 解 


首先 应 该 解 固体 外 壳 层 的 自转 角速度 Q 的 方程 式 〈18)， 然 后 再 将 Q 解 代入 


超 流 层 的 日 转角 速度 Q, 的 方程 式 C19), FERO, 的 解 。 


方程 式 〈18) 是 一 阶 线性 方程 式 ， 积 分 形式 为 


aw =e! a0 f ll - K(e# - Mie pe OO) a; 
xfe*[-(e* -p7 J pe 
令 解 的 形式 为 
QW) =e A, +1,)+Ce”. (20) 


C 为 积分 常数 。 


根据 上 面 两 式 


1 = Q(0) -fa _K(e* -1)|) “eat 
T 


可 以 证 明 积 分 式 内 第 二 项 K(e 习 一 1) <1， 故 可 以 用 二 项 式 定理 展开 ， 去 掉 第 三 项 


后 有 
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将 (21) 式 和 22) 式 代 入 (20) st 
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4 t=0， 得 到 常数 C 
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将 常数 C 代入 (23) 式 后 得 到 解析 解 
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将 上 式 O(D 代 入 〈19) 式 后 可 得 到 超 流体 内 层 角 速度 Q,。 考 虑 到 Q, est) AL 


层 的 角速度 ， 一 般 观 测 不 到 ， 而 且 推导 的 式 子 这 长 ， 故 只 研究 外 壳 层 的 自转 角 速 
度 变化 ， 略 去 了 内 层 超 流体 自转 角速度 的 解 。 


3. 脉冲 星 PSR0531+21 (Crab) 和 PSR0833-45 (Vela) 的 理论 数值 结果 


本 文 研究 具有 两 成 分 模型 的 脉冲 星 PSR0531+21 和 PSR0833-45 在 磁 豪 减 作 


用 下 外 壳 固 体 层 的 自转 角速度 变化 ， 为 此 先 列 出 两 个 脉冲 星 的 物理 参数 如 表 1。 


其 中 ou，@，zr 引 自 文 献 [5，7]， 表 面 磁场 引 自 文献 [4]，A = RB 一 般 假 定 中 


子 星 的 半径 R=1.2x10*cm， 转 动 惯量 引 自 文献 [1] T=1.4x10“ (cm - g), Wee 


MAME =1.1111x10%/yr, T =2/E=1.8x10° yr. 7! 


Sh 
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表 1 脉冲 星 PSR0531+21 (Crab) 和 PSR0833-45 (Vela) 的 物理 


Table 1 Data for PSROS31+21 and PSRO833-45 


脉冲 星 Pulsars Q, (rad/s) Q=1,/1 r(d) Bx10"(G) 


PSRO531+21(Crab) 190 0.690 7.7(d) 3.78 


PSR0833-45 (Vela) 70.5 0.145 1.2(yr) 3.38 


将 表 1 中 的 数据 和 R，LI，<& 的 数值 代入 〈11) 式 和 (13) 式 给 出 的 如 和 K 


的 数值 ， 再 取 时 间 间 隔 上 =1lyr ， 数 值 结果 如 表 2. 


表 2 脉冲 星 PSR0531+21 和 PSR0833-45 的 数值 结果 (t=1yr) 


Table 2 Numerical results for PSRO531+21 and PSR0833-45 (t=1yr) 


脉冲 星 Pulsars 4o x10 (om .GOJ) K(c-g-s) “oo anradyy rad) 


PSROS31+21 6.53 1546.2 0.9991 189.8270 -0.1710 
PSR0833-45 5.84 170.2 0.9999 70.4929 -0.0071 
4. 讨论 和 结论 


(1) 从 表 2 的 数值 结果 ， 脉 冲 星 PSR0531+21 和 PSRO833-45 在 磁 豪 减 力矩 


作用 下 壳 层 自转 角速度 逐年 随时 间 长 期 减 小 。PSR0531+21 每 年 减 小 大 于 


PSR0833-45 的 减少 ， 这 主要 是 由 两 者 的 K 值 不 同 。K 值 同 磁场 和 角速度 有 关 。 


根据 表 1, Crab 的 磁场 大 于 Vela 的 磁场 , 而 两 者 的 角速度 相差 很 大 , PSROS31421 


的 角速度 几乎 大 于 PSR0833-45 的 两 倍 以 上 , 故 前 者 的 K 值 大 于 后 者 的 K 值 , 所 


以 PSR0531+21 的 减速 大 于 PSR0833-45 的 减速 。 


(2) 同 前 文 【4】 相 比较 ， 对 于 PSR0531+21 的 磁 辐 射 作用 使 其 角 速 减 小 每 


HEIE AD = -0.2450(ra4d/) 而 在 本 文 磁 衰减 作用 使 角 速 减 慢 值 OQ. = 0.1710, 


衰减 作用 稍 小 于 磁 辐 射 作 用 。 磁 辐射 是 由 于 脉冲 星 高 速 自转 能 转换 来 的 ， 而 厂 衰 


减 是 由 于 磁场 减弱 造成 的 ,然而 自转 能 的 损失 加 快 辐射 力矩 作用 远 远 大 于 磁场 育 


减 的 作用 。 故 磁 辐 财 对 脉冲 星 自转 效应 大 于 磁 矩 衰减 对 目 转 产生 的 效应 。 


(3) 关于 (10) 式 积分 限 取 值 ， 其 下 限 t。= 0， 可 以 理解 脉冲 星 诞 生 时 的 开 


全 时间, 但 本 文 不 是 研究 脉冲 星 诞 生 时 的 演化 问题 。 另外 也 可 以 理解 脉冲 星星 震 


后 的 开始 时 间 ， 但 星 震 后 的 开始 时 间 角 动量 7, 和 Q, 值 不 好 确定 ， 而 我 们 所 知道 


的 和 @u 是 目前 开始 的 值 ， 如 表 1 的 数值 ， 所 以 在 积分 式 〈10) t 取 值 最 好 是 


从 目前 t=0 开始 为 宜 。 


(4) 关于 (21) 式 和 (22) 式 中 的 [L- 天 Ce-D] 克 和 [1 天 (e- -) 2, 


只 有 K(e“ -D<1 才 能 用 二 项 式 定 理 展开 。 首 先 估计 下 K(e”-D 的 数值 : 根据 


文 [9] Tp =1.8x10°4F E=% =1.11x10°/ 年。 本 文 取 时 间 间 隔 一 年 ，t=1 年 ， 


e% —1=0.0000011, X} F PSR0531+21,，K=1546.2，K(e -—1)=0.0017<1; 对 


于 PSR0833-45, K=170.2, K(e% -DTD=0.00018<1. 因 此 ， 两 式 可 以 用 二 项 式 定 


理 展 开 


Hke- 14 ae? ijs T jii 
Hp (e4 —1)? =1.21x1022, 所 以 第 三 项 可 以 略 去 。 


5. 脉冲 星 PSR0531+21 (Crab) 和 PSRO833-45 (Vela) 具有 星 震 的 两 个 脉冲 


星 。 前 者 每 隔 三 年 ， 后 者 每 隔 二 年 发 生 一 次 星 震 或 跃 变 。 跃 变 发 生 时 是 突然 暂 短 


加 速 后 恢复 跃 变 前 的 角速度 ， 而 这 种 跃 变 近似 周期 性 的 。 本 文 给 出 的 这 两 个 脉冲 


星 的 减速 是 长 期 角 速 减 慢 效应 ， 而 临时 突然 加 速 不 会 影响 长 期 减速 效应 。 o 


6. 本 文 得 到 的 结论 是 ， 脉 冲 星 两 成 分 模型 磁 衰 减 是 存在 的 ， 它 影响 角速度 


长 期 减 慢 的 作用 。 这 种 作用 是 通过 磁 窍 衰减 制 动 角速度 而 造成 减 慢 的 。 
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The Impact of Magnetic Decay Braking Torque on the Secular 
Retardation of Spin of Two-Components of Pulsars 


Li Linsen 
School of physics of Northeast Normal University Changchun, 130024 


Abstract: The purpose of this paper researches the action of magnetic decay on 


spin down of two-component of pulusar. The secular speed down of the magnetic 
decay braking torque on a model of two-components of pulsars is studied by using the 
analytical method. The analytical solutions of the couled equation system with the 
model of two-components are given. The analytical solutions obtained are applied to 
the research for PSRO531+21 (Crab) and PSRO833-45 (Vela). The numerical results 
show that spin of Crab and Vela pulsars speed down in 1 yr: -0.1710 (7 a. ) and 
-0.0071 (7 adi ) respectively. The results obtained are discussed in this paper and the 
conclusion that magnetic decay exists in two-component model is given especially. 
The discussions and conclusion are held on the obtained results. 
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